Themoelectric Power of High- Tc Cuprate Superconductors Calculated from the 

Electronic Structure. 
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We have calculated the thermopower of the Bi2Sr2Cu06 and Bi2Sr2CaCu2 08 superconductors 
using an ARPES-derived dispersion, with a model pseudogap, and a marginal- Fermi liquid scattering 
rate that has a minimum with respect to energy at the van Hove singularity (vHs). Good fits with 
data are achieved across the entire phase diagram, thus confirming the dispersions, the locations of 
the vHs and the dominance of the diffusion thermopower over the phonon drag contribution. 
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The thermoelectric power (TEP) of the high- 
temperature cuprate superconductors was early shown [ij 
to exhibit universal behaviour. Roughly speaking, the 
linear part of the TEP takes the form S{T) = Sq - aT. 
So is positive and large in the underdoped regime and de- 
creases with doping, becoming negative in the overdoped 
regime. This behaviour has been modelled in terms of 
separate metallic diffusion and phonon drag contribu- 
tions by TrodahlQ. In the model, the TEP comprises 
the sum of a negative linear metallic component, and a 
positive phonon drag component which rises at low tem- 
peratures before saturating at high temperatures. 

Alternatively, Mcintosh and Kaiser Q calculated the 
TEP from a model electronic structure containing a van 
Hove singularity (vHs), and were able to explain the 
TEP in terms of the diffusion component alone. Here 
the reduction in the positive part of the TEP with dop- 
ing results from the approach to, and traversal of, the 
vHs by the Fermi level {Ep). Experimental support for 
this picture has been published by Kondo et aL[J|. They 
measured both the electronic structure, by angle-resolved 
photoemission spectroscopy (ARPES), and the TEP of 
a series of Bi-2201 samples. The TEP was calculated 
from the electronic structure and a good correspondence 
with the data was obtained for over- and optimally doped 
samples. Their data and calculations (replicated by the 
present authors) are shown in Fig. [TJ The calculations 
show the positive peak in the TEP reducing and disap- 
pearing as Ep approaches and traverses the vHs. 
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FIG. 1: Measured and calculated thermopower of Bi-220l[3]- 
The curves are calculated assuming a constant mean free path. 
Dashed curves include corrections to the chemical potential 
for temperature. The dotted curve shows the effect of a pseu- 
dogap in the density of states. The values of Ef — EvHs from 
underdoped to overdoped are 70, 42, 16, 6 and -4meV. 



The diffusion thermopower is given by[5| 
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a{e) is the spectral conductivity, which under the Boltz- 
mann formalism is given by 



The calculations by Kondo et ai, which have no free 
parameters but assume a constant mean free path and 
a T-independent chemical potential (fi), are unable to 
reproduce the magnitude of the TEP in the underdoped 
regime. This is perhaps not surprising given that the 
pseudogap, a key feature of the underdoped regime, is 
not taken into account. Here we extend their approach 
by including models for the pseudogap and the scattering 
rate. 
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k 

is the mean free path given by 

£a = Wa(e,k) • T(e,k,T) 



(3) 



where Va — deih) / dka is the group velocity of the con- 
duction electrons, r(e, k, T) is the relaxation time and 
a — X or y. 
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FIG. 2: The density of states of Bi-2201 calculated from the 
ARPES-derived energy-momentum dispersion and the con- 
ductivity calculated from Eqn. [5] assuming a constant mean 
free path and constant relaxation time. 
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FIG. 3: Angle dependence of the lifetime of Bloch states at 
Ef determined by ARPES measurements at 200K[i3]. The 
OD7K data indicates long lifetimes near the saddle points at 
and 90°. 



The TEP as defined by Eq. [T] is a measure of the 
asymmetry in a{t) about /i via the spectral window 
(e — /x)[9/(e)/9e]. The saddle-points in the measured 
energy-momentum dispersion e(k)[3| give rise to a vHs 
in the density of states as shown in Fig. [21 Assuming a 
constant mean free path, a{t) oc UQ,^[e(k) — e] is just 
a velocity-weighted density of states. Because the veloc- 
ities near the saddle points are small, the peak in cr(e) is 
less prominent than the vHs is in the DOS. 

As a first step we have corrected /i so that the carrier 
concentration 



n = 2 / /(e)7V(e)de 



(4) 



remains constant with temperature resulting in the 
dashed curves shown in Fig. [T] This alTects the magni- 
tude but not the general behaviour of the thermopower. 
[1 is corrected for temperature in all subsequent calcula- 
tions. 

Next we include a model for the pseudogap that we 
have used previously to describe the electronic entropy @, 
0|, superfluid density and Raman responsejl, 01 . It is 
based on ARPES results and takes the formn 



^g^max COS 
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otherwise 



where < 6 < 7r/2is the angle subtended at (tt, tt) by 
the points {kx,ky) and (7r,0). reproduces the Fermi 
arc phenomenology and is given by[7| 



70 — VQ^max 



1 — tanh 



2T* 



(6) 



where T* — Eg^max/ks- The non-states-conserving na- 
ture of the pseudogap is implemented by removing states 
with \e{\i) - < Eg from Eqs.[Tl[2]and|31 This is equiv- 
alent to such states having very short lifetimes. 

The dotted curve in Fig. [T] shows the TEP calcu- 
lated for the most underdoped sample assuming a con- 
stant mean free path and including a pseudogap with 



E„ 



= ISmeV and ^o,; 



33°. There is a boost in 



the TEP below 250K and an increase in slope above 250K 
but the effect is too small to significantly improve the fit. 
Under the constant £ assumption the model pseudogap 
does not introduce enough additional asymmetry in a{e) 
about /i. This brings us to the scattering rate. 
We adopt a scattering rate of the form 



hT'^ (e, T) = X^iirkBTY + (e - E,hs) + a (7) 

where A is a coupling constant, a takes a fixed value 
of ImeV and is included to prevent r from becoming 
infinite. The square root term in Eq. [7] is similar to 
the implementation of the max(|w|,T) marginal Fermi 
liquid single-particle scattering rate employed by Abra- 
hams and Varma ll| , but it differs in the location of the 
minimum with respect to energy. In order to produce a 
peak in a{e) at the energy of the vHs, Eq. [7] must have 
a minimum at EyHs rather than at Ep. This is because 
^^u^(5[e(k) — e] does not exhibit a peak at E^hs, as 
shown by the plot of a{e) for constant r in Fig. [51 Sup- 
porting evidence for a small scattering rate at the vHs 
comes from angle dependent ARPES measurements of 
the lifetime of Bloch states at i^i? |10j. The data, repro- 
duced in Fig. [3l shows the average lifetime around the 
Fermi surface increasing with doping. Furthermore, for 
the most overdoped sample (filled circles) the relaxation 
time is maximal at and 90 degrees. Because Ep lies 
close to the vHs at this doping this indicates that the 
scattering rate is small near the saddle points. Similar 



3 




(a) 



100 200 300 400 

Temperature (K) 



500 




20 10 -10 -20 -30 -40 

S(293K) (|aV/K) 

FIG. 4: (a)Fits to Bi-2201 TEP data of Refs. 3 (closed sym- 
bols) and |13l ] (open symbols) using models for the pseudogap 
and scattering rate given by Eqs. (5] |6]and[7] (b) Parameters 
extracted from the fits in (a). Solid symbols are experimen- 
tally determined. 



behaviour is also found in overdoped Tl-2201[12|. The 
e — EyHs dependence arises from a more complex under- 
lying fc-dependent relaxation time and importantly in- 
troduces asymmetry in a{e) about /.t. This in turn gives 
rise to asymmetry in the pseudogap, which increases the 
TEP. 

Fits to the Bi-2201 data of Kondo et al. and Okada 
et aZ.[l3l are shown in Fig.|Ha). The parameters used in 
the fits are plotted in Fig. [4l^b) as a function of TEP at 
293K. We stress here that A is the only truly free param- 



eter used in the model. Ef 



E. 



g,7nax 



and 9o.r, 



can be determined by experiment. Experimental values 
for Ep ~ EyHs are shown by the solid blue squares. For 
the most overdoped dataset we find that a value of ImeV 
for Ep — EyHs gives a better fit than the -4meV measured 
experimentally. We note that the ARPES measurements 
were performed at 200K. The calculated chemical poten- 
tial corrections support a drift of a few meV over this 
temperature range due to the close proximity of the vHs. 

We have also performed fits to Bi-2212 TEP data us- 
ing an ARPES-derived bilayer dispersion Separate 
instances of Eq. [7] were applied to the antibonding and 
bonding bands with coupling constants Xab and Xb b re- 
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FIG. 5: (a) Fits to our own Bi-2212 TEP data (black solid 
symbols) as well as data from Refs. [H [Tsl. [Tcj . (b) Parameters 
extracted from the fits in (a). 



spectively. Recently we used the same dispersion and 
pseudogap model to fit electronic entropy data[7]. From 
those fits we were able to extract values for Ep — EyHs, 
E. 



g,max 



and 9o^max that closely matched values deter- 
mined directly from ARPES. Here we have used those 
values as a guide for fitting the TEP data, leaving only 
Xab and Xbb as truly free parameters. Fits to our own 
Bi-2212 data as well as data from Obertelli et aL[l[, 
Mandrus et al.^E^ and Munakata et al.fid\ are shown in 
Fig.mja). The corresponding parameters as a function of 
hole concentration are shown in Fig.[5l^b). 

From the fits, the curvature in the overdoped data 
above lOOK is seen to arise from the bonding band vHs, 
which is located approximately lOOmeV below the anti- 
bonding band vHs. The same conclusion was reached by 
Takeuchi et aL[l7j who performed calculations using a 
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the pseudogap on the TEP, 



6: Two of the curves from Fig. [5{a) recalculated with 
g,max Set to zero (dashed curves) to illustrate the effect of 



bilayer dispersion and a constant mean free path. 

The effect of the pseudogap on the TEP is illustrated in 
Fig.[6]for two different dopings, by setting Eg^max to zero 
while keeping Ep — E^hs and A unchanged. The pseudo- 
gap boosts the TEP at low temperatures and increases 
the high temperature slope. The peak in the TEP results 
from the combined effects of the pseudogap and the vHs, 
not just the pseudogap alone. Attempts [l8| to determine 
the pseudogap temperature T* by a simple scaling anal- 
ysis are therefore somewhat naive. The relevance of the 
band structure to the scaling behaviour of the TEP is 
also discussed by Okada et o/.[13|. 

The coupling constants are roughly constant in the 
overdoped regime and increase systematically with the 
opening of the pseudogap. The behaviour is reminiscent 
of the isotope effect coefficient [l^]- 

Recently Daou et al.^^, have measured the TEP of 
Lai.g_a;Ndo.4Sra;Cu04 and for x — 0.24 found a linear 
behavior when S/T is plotted versus log(r). Below this 
doping the low-T data (for x = 0.2) fanned upwards while 
above this doping, for x — 0.3, it fanned downwards. (In 
the latter case they used data for the Nd-free compound). 
They deduced from this the presence of a quantum crit- 
ical point Sit X — 0.24, which in their view locates the 
termination of the pseudogap line T*{x). We feel this 
conclusion is not warranted. Leaving aside the paucity of 
doping points, the behavior is just what is expected from 
Figs, m and [5] for traversing a vHs. In Fig. [7] we replot 
the latter data as S/T versus log(r). It shows a distinc- 
tive fanning out at low T either side of the vHs which 
is located between the 3rd and 4th to bottom curves. It 
should also be noted that La2-a;Sra;Cu04 exhibits a peak 
in the DOS and in the static susceptibility precisely at 
X — 0.24 indicating the location of the vHs crossing 21 1. 

In summary, we have calculated the TEP of high- Tc 
cuprate superconductors using an ARPES-derived dis- 
persion, with a model pseudogap, and a marginal-Fermi- 
liquid scattering rate that has a minimum with respect 
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FIG. 7: S/T versus log(r) for the data in Fig.EIa). Symbols 
have the same meaning. The fanning out behavior (which 
near T — cannot be scaled) is due to crossing the vHs and 
not a quantum critical point (where scaling would apply near 
T = 0). 



to energy at the vHs. Good fits with data are achieved 
across the phase diagram confirming the overall ARPES- 
derived dispersion and the location of the vHs. Our re- 
sults show that the diffusion thermopower dominates over 
the phonon drag contribution. 
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